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Anderson transition of the phonon modes is studied numer- 
ically. The critical exponent for the divergence of the localiza- 
tion length is estimated using the transfer matrix method, and 
the statistics of the modes is analyzed. The latter is shown 
to be in excellent agreement with the energy level statistics 
of the disrodered electron system belonging to the orthogonal 
universality class. 

Anderson localization, which was originally proposed 
for disordered electron system, is now widely dis- 
cussed in many random systems of classical waves, such 
as photons || and phonons. j| Especially interesting is 
the Anderson transition in three dimensional (3D) sys- 
tems, where the wave form changes from extended to ex- 
ponentially localized one. In disordered electron systems, 
this can be well described by the scaling theory, m 
and universal behaviors of localization length, ^^[con- 
ductance distribution as well as level statistics |] juj] 
have been demonstrated. The important point is that 
the Anderson transition is classified into a few universal- 
ity classes. 

Recently, Song and Kim |l3| have shown that the 3D 
phonon modes in disordered lattice undergo Anderson 
transition. However, whether this Anderson transition 
is classified into the same universality classes as elec- 
tron systems remains to be studied. In this short note, 
we demonstrate the Anderson transition of 3D phonon 
modes and discuss its critical behavior, and show that 
it is classified into the orthogonal universality class of 
disordered electron system. 

We begin with the following equation of motion for the 
lattice displacement u(r, t), 



(2) 



Ki(r)(u(r + e h t) -u(r,t))] 



where r is the lattice position, and the unit vector in 
the i-th direction (i = x,y,z). Here we have assumed 
the simple cubic structure, and the displacement is rep- 
resented by a scalar quantity for simplicity. Ki(r) is the 
spring constant connecting the site at r and that at r+e^. 
Only the nearest neighbor interaction is taken into ac- 
count. 

Putting u(r, t) — e~ iut u(r) and setting £1 = uj 2 , we re- 
duce the problem to the matrix eigenvalue problem which 
is similar to the tight binding model where the transfer 
energy as well as the potential energy is random. The 
difference is that off diagonal elements -ff^ 's and the di- 
agonal elements H^s are correlated in the phonon system 
such that 



The correlation of the matrix element is very important 
especially in the one dimensional systems. Without the 
correlation, i.e., when H^s are independent random vari- 
ables, all the states are localized. When only the off- 
diagonal elements are random, the states delocalize at 
the band center (E = 0) with the localization length 
£ ~ | log |Q. With the above correlation, the local- 
ization length diverges according to the power law of the 
mode frequency. [fL5| 

In the actual simulation, the random spring constant 
KiS are independently and uniformly distributed be- 
tween 0.1 Ko and 1.9Kq. We scale the frequency ui by 
\J K /M and the length by the lattice constant. The den- 
sity of modes D(Q,) obtained by diagonalizing 14 x 14 x 14 
systems for 935 configurations is shown in Fig. 1. 

We then discuss the critical behavior of the mode u. 
We apply the transfer matrix method (|J^] to calculate 
the decay length for N x N x L system along in- 
direction, setting L much larger than N and A at. In 
actual calculations, N is varied from 8 to 12, and the 
simulation is stopped when the relative error of A at be- 
comes smaller than 2%. 

In Fig. 2, we plot An = Xn/N as the function of f2. For 
extended modes, An is a increasing function with respect 
to N, while it is a decreasing function for localized modes. 
At the critical point An becomes N independent. The 
common crossing point in Fig. 2, therefore, gives the 
critical point f2 c , which is estimated to be 



M-^u(r,t) = J2 [K i (r-e i )(u(r-e i ,t)-u(r,t))+ (1) 



Q c = 13.0 ±0.2. 



(3) 



The critical exponent v for the localization length (or the 
correlation length) £ which characterizes the divergence 

£~|ft-O c |-* (4) 

can be estimated by expanding An as 

A N = A c + A 1 N 1 ^{n-n c ) + --- (5) 

In the present simulation, v is estimated to be 

!/»1.2±0.2, (6) 

Recent precise estimate of ^And in the Anderson model 
gives 



J/And = 1.59 ±0.03 



(7) 
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for 3D orthogonal universality class, i.e., systems with- 
out magnetic fields. Present estimate of v for phonon 
mode localization length is considerably smaller than 
;/And- This, however, should not be taken seriously, since 
when the transition occurs near the band edge, we tend 



to underestimate v. 16 

In fact, if we study the level statistics regarding f2 close 
to the critical point, we obtain the same critical level 
statistics obtained for the Anderson model belonging to 
the orthogonal universality class. The distribution func- 
tion P(s) for the nearest neighbor spacing of the two suc- 
cessive modes s near the critical points (12.5 < ft < 13.5) 
is plotted in Fig. 3, together with the critical level statis- 
tics for the orthogonal universality class in electron sys- 
tems. H The coincidence between the two is excellent, 
and we conclude that the Anderson transition of phonon 
modes is classified into the orthogonal universality class 
in the disordered electron systems. The critical behav- 
ior in 3D is not altered by the correlation of the matrix 
elements (Eq.(||)). 

In the present model, we have assumed the simple cu- 
bic lattice where the general vector displacement prob- 
lem can be reduced to the scalar displacement problem 
(isotropic Born model). |l7j In real lattices, all the com- 
ponents of the displacement vectors should be taken into 
account. @@ Such a model is very difficult to analyze 
by the transfer matrix calculation. The present method 
of level statistics might be a powerful tool to investigate 
the phonon Anderson transition in such realistic lattices, 
which is a problem left for the future. 
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Figure captions 

Fig. 1: The density D(Vl) of the square frequency f2. 

Fig. 2: Normalized decay length A^r vs.O for N =8 
(A), 9(o), 10 (□) and 12 (o). The dashed lines are the 3rd 
order polynomial fit to the data. Common crossing point 
indicates the critical point f2 c , and the A-dependence of 
the slope at VL C gives the localization length exponent v. 

Fig. 3: Distribution function P(s) for the nearest 
neighboring STs close to f2 c (•). That for the disordered 
electron systems (o) as well as the Poisson (broken line) 
and the Winger (solid line) distributions is also shown 
for comparison. 
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